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channel, is studied. It is found that those models have either Sp(2s + l)-symmetry for the fermions 
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INTRODUCTION 

Recently, the study on cold atoms with high spin has 
aroused lots of attention in the fields of atomic, molec- 
ular, optical and condensed matter physics. Due to the 
spin exchange interactions, many interesting spin ordered 
states arise and the phase diagrams of these systems 
are very rich. For instance, in the spin-1 spinor Bose- 
Einstein condensations, the bosons are found to form 
the pairs and the pairs condense even in the repulsive 
regime [iHl]- In experiments, by using the atom cool- 
ing and trapping techniques, one can prepare the high 
spin cold atomic systems, such as ''Li ^"^Na ®^Rb 
0-0] with hyperfine spin 1; ^'^Cr [§] with hyperfine spin 
3/2; and 0, i^syb [uj, ^Ca [H, »^Sr [H, i^^Cs 



with more higher ones. Using the Feshbach res- 
Til . 17 1 and confinement induced resonance 18 1 



onance 

technique, the interactions among the atoms can be ma- 
nipulated. This provides a good platform for studying 
the tunable condensed matter systems. In theoretical 
approaches, the low-energy effective models of the dilute 
ultracold atomic systems are the quantum gas with con- 
tact interactions, and the spin exchanging interactions 
should also be considered for systems with internal de- 
grees of freedom [3, ISl . 

Symmetric analysis plays a very important role in 
studying the quantum many-body systems. Physical 
properties such as ground state manifold and order pa- 
rameters are closely related to the symmetry of a system 



21|-|23j . The analysis of the symmetry can give us some 
hints to do the suitable approximation and study the 
physics such as phase diagram in the frame of mean-field 
theory. It can also simplify the analytical and numeri- 
cal calculations. In the cold atomic systems with delta 
function interactions, the symmetric or anti-symmetric 
properties of the identical particles restrict the forms of 
spin exchange interactions. Effective spin exchanging in- 
teractions only take place in the channels with symmetric 
spatial wave function. Such property may make the sys- 



tems to have intrinsic symmetries in the spin sector. For 
example, in the spin-3/2 system, the S0{5) symmetry is 
found 

The strong quantum fluctuation and correlation make 
the physics of one-dimensional (ID) system quite differ- 
ent from the ones of higher dimensions. Many numerical 
and analytical methods are developed to study the ID 
systems. The exact solution is a good starting point to 
study these systems, since it can give us conclusive re- 
sults. A well-known exactly solvable system in the cold 
atomic system is the Lieb-Liniger model 2S], where the 
scalar bosons are studied. The fermion case with spin-1/2 
are exactly solved by Yang [2^. Sutherland generalized 
Yang's model to arbitrary spin case, where the system 
has the SU (2s -I- 1) symmetry [i^], with s the spin of the 
particles. At present, it is clear that the multicomponent 
quantum gas including the Bose-Fermi mixtures with 
delta function interactions are integrable, if the masses of 
each species are equal and the interactions are equal [28- 
[30j |. In these models, the scattering lengthes in different 
channels are the same, for that the spin exchange inter- 
actions are not considered. However, the spin exchanging 
usually can not be neglected in experiments, and many 
novel ordered states are induced by the spin exchanging. 
Motivated by this consideration, we proposed a SO {3) 
integrable spin-1 bosonic model 13 ll and a 5^(4) inte- 
grable spin-3/2 fermionic model [32|, where the contact 
spin exchange interactions are considered. The scattering 
lengthes in different channels are different. 

This paper is a generalization of our previous works 
[31] , 32 1 to cold atom systems with arbitrary hyperfine 
spins. For a special interaction form of atoms with hy- 
perfine spin s, we find that the fermion system has the 
U{1) (E) Sp{2s + 1) symmetry while the boson system has 
the U{1) (Xi 5*0(25 -I- 1) symmetry. The generators of 
corresponding algebra are constructed by the magnetic 
multipole operators. Based on the symmetry analysis, 
we propose a new class of exactly solvable models in one 
dimension. 
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(a) scattering process of P„i (b) scattering process of 



FIG. 1: The sketches of two-atom scattering processes. 



MODEL AND SYMMETRY 

For the delta function interaction models of dilute cold 
atomic gas with hyperfine spin s, the spin exchange in- 
teraction between two particles i and j is usually writ- 
ten as spin projection operators P/^- in different channels 
with total spin-/ (/ = 0, 1, 2, • • • , 2s). Nontrivial scatter- 
ing processes occur only in the even I channels because 
of the symmetry or anti-symmetry behavior of the wave 
functions. To study the behavior of such systems away 
from the SU (2s-f 1) symmetry point, we consider a simple 
case, i.e., all the scattering lengths of nonzero I channels 
are the same. The Hamiltonian reads 



N 



H 



N 

-E 



'C2 J2 

1=2,4,- 



PL 



d{n~r,). (1) 



Here, N is the number of atoms, is the position of the 
i-th atom, ci is the interaction strength in spin-0 channel 
and C2 is the one in other channels. 

The two-body scattering in the system ([TJ is quite in- 
teresting. There are two kinds of scattering processes 
in the spin sector. One is |s, m; s, m')(s, m'; s, m| + h.c. 
provided by two-particle permutation as shown in Fig. 
[Ua), here m and m' are the spins along z direction, 
and m,m' — s, s — 1, • • • , —s. In this process, the par- 
ticle numbers with different m are invariant. As a 
consequence, the total spin S and total particle number 
N are also conserved. The other scattering process is 
|s, m; s, —m) {s, m'; s, — m'| -|- h.c. provided by the projec- 
tor operator P^^^, as shown in Fig. [Ijb). In this pro- 
cess, two particles with opposite spins scatter into an- 
other pair, and the absolute value \m'\ can be unequal 
to |m|. Obviously, this process does not affect the total 
particle number iV, but it destroys the invariability of 
Nm- Therefore, the particle number N„i no longer con- 
serves. Nevertheless, by careful consideration, we find 
that Jm = Nm — N-m with m = s, s — 1, • • ■ and m > 
are still invariant. The invariance of Jm means that the 
total spin S = J2m ''^Jrn is still a good quantum num- 
ber. Besides, some of the magnetic multipole operators 



are also invariant. The multipole operators are observ- 
able physical quantities and can be defined in the form 
of irreducible tensors. 
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fl - ^{l + m)\{l^m)\/{2l)\(J{ J2 s'™Jd:t; ™, ' 

j=l nij 

m = -I, -I + !,■■■ J, / = 1,2,- • • ,2s, 

l_i = (r-is^)/V2, lo = V2s^ li = -(s"+is^)/V2. 

(2) 

Here, s" {a — x, y, z) is the spin operators of one particle, 
rrij = —1, 0, 1, and the sum E' means I'm — '^^ mi\ < j 
for any j. The total multipole operators for N particles 
are iW^ = Si^i ^m'; where T/jI' is the /-rank multipole 
operator of the i-th particle. It can be proved that the 
multipole operators with odd rank are commutative with 
the Hamiltonian 



0, 



(3) 



where V = I]^j (ciP°- -I- C2j2i=2,4,-- - Hj)^ and thus are 
the conserved quantities of system ([1]). This can be un- 
derstood from the two-body scattering processes. If we 
only consider the process of permutation, all multipole 
operators are commutative with the spin part of the 
Hamiltonian, since the exchanges of spins have no ef- 
fect on the magnetic properties. However, the scattering 
processes of make the magnetic quadrupole change. 

The odd rank multipole operators can be used to con- 
struct the generators of Sp{2s + 1) (half odd s) and 
S'0(2s + 1) (integer s) algebras. Because the algebras 
so(s, s + 1) and so(2s + 1) have the same complex ex- 
tensions so(2s -I- 1,C), and if a model possesses the 
SO{s, s -f 1) symmetry, it must have the S0{2s + 1) 
symmetry. Here we use group SO{s, s + 1) instead of 
S'0(2s -I- 1) to reveal the hidden symmetry in the system 
(dl). The generators of these groups can be defined uni- 
formly as the (2s + 1) x (2s + 1) matrices satisfying the 
following conditions 



try 0, r 



-Y* + jy = 0. 



(4) 



Although the operator J is different for Sp{2s + 1) and 
SO{s, s -f 1) groups, we can write it in a uniform repre- 
sentation 



7" 



i-iy 



m, m 



s, s ■ 



(5) 



Note that it is not valid for SO{s, s -f 1) group with half 
odd s, and we use SO{s, s + 1) group for the ones with 
integer s (i3„ type algebra, n — s) in the following dis- 
cussions. 

Eq. (m shows a similarity of Sp{2s+1) and SO{s, s+l) 
groups. This enables the different symmetries of Hamil- 
tonian ID) to be proven simultaneously. Since the mul- 
tipole operators defined in eq. are all (2s + 1)^ — 1 
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linearly independent (2s + 1) x (2s + 1) real matrices with 
zero trace, we can use them to construct the generators 
of Sp{2s + 1) and SO{s, s + 1) group. If we choose the 
basis of multipole operators ^ as the eigenstates of s^ 
and the representation of J takes the form of (O, the 
generators of Sp{2s + 1) and SO{s, s + l) groups defined 
in (|4]) can be represented as 



eters cq and C2' 

AT „, N 



(7) 



Yi = ITl, (odd /). 



(6) 



Therefore, (s + l)(2s + l) operators Yj^ are the generators 
of Sp{2s + 1) group if s is half odd, and the s(2s + 1) 
operators Y^ are the generators of SO{s, s + l) group if 
s is integer. For that multipole operators with odd rank-^ 
are commutative with the spin part of Hamiltonian, eq. 
([3|), the corresponding symmetries hold for the system 
([!]). There is also a J7(l) symmetry for the coordinate 
part, then the system ^ has ?7(l)(8)>5'p(2s + l) symmetry 
for the fermionic case and U{1) (E) SO{2s + 1) symmetry 
for the bosonic case. 

There are three homomorphisms Sp{2) ~ SU{2), 
SO{3) ~ SU{2) and Sp{A) ~ SO{5). For the case 
s — 1/2, only one channel P^j^ involved, the model in 
the spin sector has the SU{2) {Sp{2)) symmetry. For the 
cases s = 1 and 3/2, two channels P°f, and P^^ are in- 
volved. When s 1, the model has SU(2) (50(3)) sym- 
metry, and when s = 3/2 the system has SO{5) {Sp{4)) 
symmetry, which are consistent with the results obtained 
in the references 
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When ci = C2, the symmetry of the system ([T} in the 
spin sector degenerates into the SU{2s + 1) one, where 
all the interaction strengthes in different channels are the 
same. The interaction of the spin part is the spin permu- 
tation operator up to a constant. The permutation op- 
erator acting on the symmetric wave functions gives the 
eigenvalue 1, and the one acting on the anti-symmetric 
wave functions gives — 1 . Thus the effective interaction is 
just the contact interaction and all magnetic multipoles 
are conserved. This can also be explained from the view 
that only the permutation operators are involved in the 
scattering process. In the form of multipole operators, 
the (2s + 1)^ — 1 generators of SU{2s + 1) group read 
Yf„,^i"'ifl+flJ/2 and Y,'J„ =i™+i(f^-fi„)/2. 



With standard coordinate Bethe ansatz method, the 
wave function of the system ^ is assumed as 



QV 



M'P)e 



<.8) 



Here, is the spin component along z-direction of i-th 
particle, nii — s, s — 1, • ■ ■ , — s, and ki (i — 1,2, ■■■ , N) 
are the quasi-momenta carried by the particles. Q and 
V are all A^! permutations of {1,2,- •• ,N}, and Qi 
( Vi) is the i-th number of the permutation Q (V). 

^{^Qi < ■ ■ ■ < a^Qw) — I\i=2^i^Q^ ~ ^Qi-i) is continu- 
ous multiplication of step function 9{x). When x > 0, 
9{x) = 1, and otherwise 9{x) = 0. Thus the function Q 
divides the coordinate space into iV! intervals. 

The two-particle scattering occurs at the interface of 
two adjacent coordinate intervals Q and Q' . The scatter- 
ing matrix of particles a and b carrying different quasi- 
momenta is defined in the two-particle spin space to de- 
scribe the relation of the superposition coefficients 

A{Q,V) = Sa,^a,^,{k^,^-k^,)A{Q',V'), (9) 
Q' = {Qi, Q2, • ■ ■ , Qc-i, Q«+i, Qe, Qc+2 • • • , Qjv}, 
V = {Vi,V2r-- ,V^-i,Vi+i,V^,V^+2--- ,Vn}, 

where ( = 1, 2, • • • , TV - 1, = a, Q^+i = b and 
A is the vector denotation of superposition coefficients 
Ami.--- ,mN ■ In the system ([7]), the wave function should 
be continuous and the first-order derivative of the wave 
function with respect to coordinate should be discontin- 
uous. Solving the Schodinger equation and using the 
symmetry or antisymmetry condition, we can obtain the 
scatting matrix. For the Sp(2s + l)-invariant fermionic 
model, the two-body scattering matrix is 



2s 



1=1,3, 



^ A 

l=2,A.--- 



-IC 



A-(s + §)ic ^o 
A-f(s + |)ic 



EXACTLY SOLVABLE MODELS 



For the iS'0(2s -I- l)-invariant bosonic model, the scatter- 
ing matrix reads 



In one dimension, it is well-known that at the iS'J7(2s-|- 
1) symmetry point, the model is integrable. As we 
showed in the spin-1 [sij and spin-3/2 cases, there 
is indeed another integrable point. For the Sp{2s + 1) or 
S'0(2s + l)-invariant Hamiltonian ([IJ, we construct the 
following exactly solve model by constricting the param- 



V- p; , V- ^-icpi A+(s+ l)ic 



The scattering matrices ([T0| and (fTTj) are different. 
In order to prove the integrability of the bosonic and 
fermionic models uniformly, we introduce the i?-matrix 
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for these two kinds of symmetries by the fohowing map- 
ping, 



'^'^ ' \[a{X)-b{Xm,sil\~X), (integers). ^ ' 
With this mapping, the explicit form of i?-matrix is 



(half odd s). 



R^:,\X) = biX)I + a(A)P,, - (2s + 1)z''>{X)p::,. (13) 

Here, 6(A) — ic/(A + ic), a(A) — A/(A + ic), / is a unitary 
operator and z''*'(A) is a scalar function depending on s, 

z(^)(A) = (-l)2^5(A)a(A/[,s + 1/2 + {-1^+']). (14) 

After some calculations, we find that for any spin-s, inte- 
ger or half odd, i?^*-* satisfies the Yang-Baxter equation 

Rab{X)Rbc{X + l^)Rab{y) - Rbcil^RabiX + m)4c(A)(15) 

In the derivation, the following relations have been used 

2s 

P°,n, = (2,s + 1)P°,P°„ Pab = Y^{~lf'-'PU. (16) 

;=o 

Since there are two invariant mappings, i?(A) i— > 
f{X)R{X) and R[X) i-> i?(— A), for the Yang-Baxter equa- 
tion [sst . Hamiltonian (jT]) is integrable. 

i?-matrix defined in eq. (fT3|) only has two sets of so- 
lutions of the Yang-Baxter equation (jisp . One set is 
z'-*-' — where the system has the SU{2s + 1) symme- 
try. In this case, there are no effective spin exchange 
interactions. The other set is eq. (|14p . The system 
has the Sp{2s + 1) symmetry for half odd s and the 
SO {2s + 1) symmetry for integer s. The corresponding 
integrable spin chains are Kennedy-Batchelor models in 
(34| . When s = 1/2, the S'p(2)-invariant integrable model 
is discussed in [26|; when s = 1 , the S'0(3)-invariant inte- 
grable model is discussed in [3l| ; and when s = 3/^ the 
S'p(4)-symmetry integrable model is discussed in [32|. 

The integrable model ([7]) has one tunable interacting 
parameter c. For the Sp{2s + 1) fermionic model, the 
interaction is repulsive when c > and is attractive when 
c < 0. For the SO{2s + l) bosonic model, the interaction 
in spin-0 channel is attractive and that in other channels 
is repulsive when c > 0, while the interaction in the spin- 
channel is repulsive and is attractive in other channels 
when c < 0. To obtain the exact energy spectrum of 
the system, we need to determine all the values of quasi- 
momenta kj . This can be done by solving the eigenvalue 
problem given by the periodic boundary condition, in 
which we can obtain the Bethe ansatz equations. 



EXACT SOLUTIONS 



algebraic Bethe ansatz method. The Bethe ansatz equa- 
tions of integrable quantum gas models are composed of 
the ones of the coordinate part, i. e. U{1) symmetry and 
the ones given by the spin part. The spin sector usually 
has nesting integrable symmetries for high spin models. 
Using the method suggested in [3^ [s^, we can obtain 
the Bethe ansatz equations for the Sp{2s + 1) model (C„ 
type algebra, n = s + 1) and the SO{2s + 1) model. 

For 5*^(25-1- 1) (s > 1/2) case, there are s-l-3/2 sets of 
coupled equations. When s > 3/2, the equations are 



n 



=1 "-J 



n 



1 Aj ^\ -12 



n 



^'-Af^'- 
1 Aj ' 



j = l,2,... ,iV, (17) 



-1^, ^^'A<'>-A<','+ic 



n- 



'2 



Z = s,s-1,... ,5/2, j = l,2,... ,M('), (18) 



■Q Aj ~\ 



1 ^f-Af -i 



• c ^ (f) ^ (i) 

n ~ 



ic ^Af'-Af+ic 



A1 



n 



A7 -ic j,_,^A/-A/,-ic 



■Xt -Xf H-ic 



»'"At-A*+2k 



i=i AJ - A/ -ic 



n 



Ay-A|'-2ic 



j = l,2,.-. ,M('), (19) 
, j = l,2,... ,M(').(20) 



Here, M^) is the numbers of rapidity A^'^, M^^^+i) =Ar^ 
and A^"^^'' = kj. When s = 3/2, the Bethe ansatz equa- 
tions degenerate into the ones obtained in [s^l- When 
s = 1/2, the system ([7]) degenerates into the Sp{2)- 
invariant spin- 1/2 Fermi gas, and the Bethe ansatz equa- 
tions are given in (26| . 

For the SO{2s + l) bosons, the Bethe ansatz equations 
have s -\- 1 sets, and when s > 1 they are 



lii.,_K_ic-l--l- 



Af(i + i)^(i)_^(! + i)_; c Ai('-i),^(i)_ 



n 



-1^ 



n 



^15 a;"- A' 



((-1) 



j = 1,2, •••,7V, (21) 

AfW 1,(0 • 

=n ~ 



;^A--A-+ic^ 



s, s - 



A/(2) •>(!) >(2) 



n 



A/(i)^(i)_^(i) 



..r^^'-^r+- /iA^>-A-+ic' 



,2, j = l,2,' 



1,2, 



,M<'),(22) 
,M('). (23) 



When s = 1, the above Bethe ansatz equations degener- 
ate into ones obtained in 31 1 . 

Therefore, if the quasi-momenta fc's satisfy the Bethe 
ansatz equations, ® is the eigen- wave- function of the 
system and the corresponding eigenvalues of energy and 
momentum are 



For the integrable systems with high symmetry, the 
exact solutions are usually obtained by using the nested 



N 



N 



E : 



K 



(24) 
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The total spin is S' = siV - ^l; M^^^ ■ 

Obviously, the Bethe ansatz equations of the present 
system are different from the SU{2s + 1) ones. The 
physical properties can be obtained from the solutions 
of Bethe ansatz equations. For example, solutions of 
50(3)-invariant spin-1 bosonic model show that there 
are bound states in the regimes of c > and c < [si} . 
for that there always exist attractive interactions in some 
scattering channels. 



REPULSIVE FERMIONS 

For the repulsive fermionic models, detailed analysis of 
the Bethe ansatz equations shows that all quasi-momenta 
k are real, which means there are no charge bound states, 
and the spin rapidities A^''' form strings. In the thermo- 
dynamic limit, the string solutions read [35| 



I 



a) spin-5/2 





FIG. 2: The ground state dressed energies of spin-5/2 (a) and 
spin-9/2 (b). Here c = 1 and n = 1. 



Ai'l,, = +(„ + !- 2j)ic/2, J = 1, 2, • • • 



3/2 </ < s, (25) 

Ai^S = ^nf + (n + 1 - 2j)ic, J = 1, 2, • • • , n. (26) 

Here, Xn]z denote the real parts of the n-string rapidities, 
and M,!'^' is the number of n-strings 



1,2, 



for A*-'-*. Based on the above string hypothesis, the finite 
temperature thermodynamic properties of the system can 
be obtained. If the temperature tends to zero, only the 
real rapidities and 2-strings for A'-'' (3/2 < Z < s) are left 
in the ground state. Substituting these solutions into the 
Bethe ansatz equations and taking the thermodynamic 
limit, we obtain the coupled integral equations. Solving 
these equations, we obtain the numbers of the z-string 
A*^'-' analytically 



M- 



l- 



2(,s 



^±1-N,{l>^-),M^^J-N. (27) 



Thus the numbers of A^'' are M(') = N{1 > 3/2), Af(^) = 
N/2, and the conserved quantities Jm = in the ground 
state. The total spin is zero, so that the ground state 
is spin singlet state. Since the string distributions are 
symmetric around the real axis, the total momentum K 
of the grounds state is zero. 

The dressed energy of charge rapidities k in the ground 
state satisfies the following equations. 



e(fc) = fc2 



^(^+i)*e(fc), |fc|<Q, 



{k)^DfKe{k), 



2' 



(28) 



Here fi is the chemical potential, * is an integral operation 
defined by * y{x) = J w{x — x')y{x)dx' , Q is the Fermi 
point which is determined by the particle density n — 



N/L, and the kernels D{k) of integral operators D are 



D{k) = ai/2 * D['\k) + ai * D\'>{k), 



is), 



1 



sin( 



2l~l , 
2s+l ' 



2s + 1 



cosh( 



27rfc/c 
2S-I-1 



cos( 



21-1 , 



(29) 



(s-i+3/2) 



* D 



(1/2) 



(fc), 3/2 </ < s. 



i:»(i/2)(fc) ^ i/[(25 _^ 3) cosh[7rfc/(c(s + 3/2))]], 

where a{x) = i/[7r(a:^ -f [tcf')]. The dressed energy for 
c — 1 and n = 1 is shown in Fig. [2] 

The physical properties of such ID systems are con- 
trolled by the parameter 7 = c/n [25[. When 7 — cx), 
we obtain the density of states, energy and Fermi point 
in the strong repulsive limit as 

p{k) = ^,{k<\Q\), E^^Q\ Q^nn. (30) 

When 7 — )• 0, the system degenerates into the free 
fermions and we have 



27r 



37r 



CONCLUSION 



2s + l 



(31) 



In conclusion, we find that there is a hidden symme- 
try of the high spin cold atomic systems with a special 
interaction form away from the SU{2s + 1) symmetry 
point. Based on the symmetry analysis, a new class of 
integrable models for cold atoms with arbitrary spin is 
proposed. 
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